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INFLUENCE OF DISCRETE RING STIFFENERS AND PREBUCKLING 


DEFORMATIONS ON THE BUCKLING OF ECCENTRICALLY 
STIFFENED ORTHOTROPIC CYLINDERS* 

By David L. Block 
Langley Research Center 

SUMMARY 

The buckling of eccentrically stiffened perfect orthotropic cylinders is the subject 
of an analytical investigation in which prebuckling deformations and the influence of dis- 
cretely located ring stiffeners are taken into account. By using energy principles, non- 
linear prebuckling equations and linear buckling equations are derived. Solutions for 
simply supported and clamped boundary conditions and for loadings of axial compression 
and internal or external pressure are obtained by the method of finite differences and a 
modified Gaussian elimination technique. For comparative purposes solutions are also 
presented for the case where the ring properties are averaged and a Galerkin solution is 
presented for the classical buckling (constant prebuckling deformations) of a simply sup- 
ported cylinder with discrete rings. Computed results for two types of simply supported 
large-diameter stiffened cylinders loaded in axial compression illustrate the influence on 
the buckling load of prebuckling deformations, discrete rings, and eccentricity of applied 
axial compressive loads; results for clamped longitudinally stiffened cylinders are com- 
pared with previous test results. 

The results show that the predicted buckling load may be substantially affected by 
consideration of prebuckling deformations, eccentric loading, and discreteness of rings 
and that agreement with experiment may be improved by including these factors. 

INTRODUCTION 

In the design of lightweight structures, an understanding of the various effects of 
stiffening elements on the buckling behavior of cylindrical shells is important. The 
importance of one of these effects, the eccentricity or one-sidedness of stiffening ele- 
ments, on the buckling strength has been demonstrated in several studies. (See 

*The information presented herein is based upon a thesis offered in partial fulfill- 
ment of the requirements for the degree of Doctor of Philosophy in Engineering 
Mechanics, Virginia Polytechnic Institute, Blacksburg, Virginia, June 1966. 



I 


I I I I I II 


refs. 1 to 8.) However, differences between experiment and theory (ref. 7) show a need 
for additional considerations in the theories, such as prebuckling deformations, discrete 
stiffeners, initial imperfections, and boundary conditions. 

The importance of prebuckling deformations has been demonstrated for unstiffened 
isotropic cylinders in references 9 and 10 and for stiffened isotropic cylinders in refer- 
ence 6. In reference 6, however, only simple supports are considered and an approxi- 
mate one-mode Galerkin solution is used to obtain the buckling loads. The importance 
of discrete rings has been studied in reference 11 for simply supported orthotropic 
cylinders in axial compression, but the theory used did not consider eccentricities of 
stiffeners or prebuckling deformations. In reference 5 a discrete-ring theory is pre- 
sented for stiffened isotropic cylinders but prebuckling deformations are neglected. 

The purpose of this paper is to present a theoretical investigation of the effects of 
both prebuckling deformations and stiffener discreteness on buckling of eccentrically 
stiffened cylinders. The effect of discrete stiffeners is studied only for rings because 
in stiffened cylinders of practical size, the longitudinal stiffeners or stringers are 
closely spaced and therefore can be assumed to be averaged or "smeared out" over the 
stringer spacing. The cylinders studied are assumed to have simply supported or 
clamped boundary conditions and to be loaded by axial compression and a constant pres- 
sure. Also investigated is the influence of eccentricity of the applied axial loads; that is, 
the effect of applying the axial-load resultant at a surface other than the cylinder's mid- 
dle surface. (See, for example, ref. 6.) 

The analysis presented is applicable to a broad class of cylindrical shell struc- 
tures, because the shell considered is constructed from orthotropic material with stiff- 
ening elements on its surface. By proper definition of the orthotropic material constants 
such an analysis can be used to predict buckling for a wide variety of stiffened shells - 
for example, sandwich-type, corrugated, filament-wound, or isotropic cylinders. 

Nonlinear prebuckling equations and boundary conditions and linear buckling equa- 
tions and boundary conditions are derived for cylinders with discrete rings (appendix A). 
Solutions to the prebuckling and buckling equations are obtained for simply supported and 
clamped boundary conditions by the method of finite differences and a modified Gaussian 
elimination procedure. (See, for example, refs. 12 and 13.) For comparison, solutions 
are also presented for cylinders with rings that are considered to be smeared out (appen- 
dix B), and for the classical buckling problem of a simply supported cylinder with dis- 
crete rings and a membrane stress state prior to buckling (appendix C). Because of the 
large number of parameters involved in the problem, results of a general nature are 
impractical to present; thus, results in the form of sample calculations are presented for 
two large- diameter axially compressed simply supported cylinders, and for clamped 
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longitudinally stiffened cylinders of the same proportions as the test cylinders of 
reference 7. 


SYMBOLS 


ApBjjCj 


A,B,C 


D x’ D y 


cross-sectional area of stiffener 


4x4 matrices defined in appendix D 


bending stiffness of orthotropic plate in longitudinal and circumferential 
directions, respectively 


xy 


E 


E ,E 
x’y 


E. 


F 1’ F 2> 


G 



I 


J 

K,L 


longitudinal bending- stiffness parameter, 


D. 


E l E A z_ 2 
s s , s s s 


1 - V, 


twisting stiffness of orthotropic plate 


Young’s modulus 


extensional stiffness of orthotropic plate in longitudinal and circumfer- 
ential directions, respectively (for corrugated cylinders see ref. 8) 


longitudinal extensional-stiffness parameter, 


S E A 

_x s s 


1 " y 


.,Fjj constants defined in appendix D 
shear modulus 

in-plane shear stiffness of orthotropic plate 

moment of inertia of stiffener about its centroid 

moment of inertia of stiffener about middle surface of shell 

torsional constant for stiffener 

matrices defined by equations (28) 
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Oil 


Li 

M 


M x>My> M xy 

M x> M y> M xy 

N x’ N y’ N xy 


N N N 
x’ x y’xy 


N„ 


N 


P i’Qi 

R 

S 


u,v,w 


matrix defined by equation (27b) 

buckling- moment function of x defined by equation (35) 
middle-surface moment resultants for stiffened shell 
middle-surface moment resultants for orthotropic shell 
middle- surface stress resultants for stiffened shell 
middle-surface stress resultants for orthotropic shell 
magnitude of externally applied axial compressive load 
total number of rings on cylinder 

4x4 recurrence-relation matrices defined by equations (44) and (45) 
recurrence- relation matrices defined by equations (31) and (33) 
radius of cylinder to middle surface of orthotropic shell (see fig. 1) 
matrix defined by equation (27c) 
matrices defined by equation (28) 

buckling-displacement functions of x defined by equations (35) 



r u 



M 




a length of stiffened cylinder (see fig. 1) 

a a’ b n’ c n series coefficients of equations (C4) 

4 4 4 
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series coefficients of equations (C6) 


a .b ,c 
r’ r’ r 


e 

h 

j 

k 

l 


m 

n 


P 

<i,r,s 

t 

t 

u,v,w 


x,y,z 


z 


a v a 2 

a 


stringer spacing (see fig. 1) 

distance from middle surface of cylinder to line on which N x acts 

one-half the height of the stringer section (see fig. 2(b)) 

integer 

total number of finite-difference stations along length of cylinder 
ring spacing (see fig. 1) 

number of finite-difference spaces A between adjacent rings 

number of full waves in cylinder buckle pattern in circumferential 
direction 

pressure (positive for external pressure) 
integers 

thickness of isotropic- cylinder shell wall 

A s 

effective wall thickness of stiffened isotropic cylinder, — + t 

d 

middle- surface displacements of orthotropic shell in x-, y-, and 
z-directions, respectively 

orthogonal curvilinear coordinates with origin lying in middle surface of 
orthotropic shell (see fig. 1) 

distance from centroid of stiffener to middle surface of orthotropic shell 
(see fig. 1), positive if stiffener lies on external surface of shell 

constants defined by equations (B3) 

4x4 selection matrix defined in appendix D 
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r l’ r 2 ,r 3’ r 4 

r xy 

r xy T 

A 

6(x-jZ) 


e x’ e y 

e x T ,e y T 

e y 

/f 

ApAg,. . v -^24 

M x ,My 


M x ,My 


constants defined by equations (9) and (23) 
constants defined by equations (B3), (B4), and (B5) 
shearing strain at middle surface of orthotropic shell 
shearing strain in orthotropic shell 
distance between adjacent finite- difference stations 

p+OO 

Dirac delta function defined so that l f(x)<5(x-jZ)dx = f(j<0, where 
6(x-jZ) = 0 when x # j£ 

Kronecker delta, equal to 0 when i * j and equal to 1 when i = j 
normal strains at middle surface of orthotropic shell 
normal strains in orthotropic shell 
strain in rings and stringers, respectively 
constant defined by equation (B3) 
constants defined in appendix D 

Poisson's ratios for bending of orthotropic plate in longitudinal and cir- 
cumferential directions, respectively 

Poisson's ratios for extension of orthotropic plate in longitudinal and 
circumferential directions, respectively 
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II total potential energy of stiffened cylinder 

n c strain energy of orthotropic shell 

n L potential energy of external forces 
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n r> n s 


$ $ 
r’rq 

v& \Er 

r’^rq 


Subscripts: 


strain energy of rings and stringers, respectively 


trigonometric series defined by equations (C8) 


A 

B 

H 


prebuckling state 

small changes away from prebuckling state which occur at buckling 

integer denoting finite-difference station at which ring is located (see 
eq. (19)) 


i 


index of finite-difference station 


max 


maximum 


r rings (circumferential stiffening) 

s stringers (longitudinal stiffening) 


A subscript preceded by a comma indicates partial differentiation with respect to 
the subscript. 


Primes indicate total derivatives with respect to x. 


GOVERNING EQUATIONS AND BOUNDARY CONDITIONS 


To develop the buckling theory for the stiffened cylinder, shown in figure 1, several 
basic assumptions are made. The stiffened cylinder is composed of an orthotropic shell 
stiffened by uniform, equally spaced rings and stringers, all having elastic properties. 
The elastic constants of the orthotropic shell are taken as those given in reference 14, 
with the assumption that the transverse shearing stiffnesses of the shell are infinitely 
large. The stringers are assumed to be closely spaced so that their elastic properties 
can be averaged over the stringer spacing. However, the rings are considered to be 
located discretely along circumferential lines on the cylinder. In cases where rings and 
stringers lie on the same surface of the shell, the effect of the joints in the stiffener 
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framework is ignored. The only applied loadings considered are an axial compressive 
end load and a constant pressure load. 

The equations are derived by obtaining strain-energy expressions for the shell, 
stiffeners, and loadings and applying the method of minimum potential energy to obtain 
the nonlinear equilibrium equations and boundary conditions of the configuration. This 
derivation is presented in appendix A and is similar to the derivation presented in refer- 
ence 4 except that orthotropic constitutive relations are employed for the shell and the 
rings are treated as discrete members. 

The nonlinear equilibrium equations which govern an eccentrically stiffened ortho- 
tropic cylinder with discrete rings are 


N v 

+ N„. = 0 

x,x 

xy,y 

N „ 

+ N „ „ = 0 

y>y 

xy,x 

-M 

+ M 

x,xx xy,xy 


M y,yy + r ' N x w ’xx ' V’yy ‘ 2N xy w -xy + *> ' 0 


(i) 


where a comma denotes partial differentiation with respect to the subscript. The bound- 
ary conditions to be satisfied at each end of the cylinder are 


N x + N x = 0 


N xy = ° 


M x + N x® = 0 


or u = 0 

or v = 0 

or w, = 0 


M x,x - M xy,y + N x w ’x + N xy w -y = 0 or w = 0 


( 2 ) 


J 


The stress and moment resultants in equations (1) and (2) are given by equations (A14) 
in appendix A. 

The nonlinear equilibrium equations and boundary conditions (eqs. (1) and (2)) are 
used to obtain the equations and boundary conditions which govern the prebuckling and 
buckling of a stiffened cylinder with discrete rings. To obtain the prebuckling and buck- 
ling equations, the displacements of the shell u, v, and w are assumed to be separa- 
ble into two parts as follows: 
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u(x,y) = u A (x) + u B (x,y) 
v(x,y) = v A (x) + v B (x,y) V 
w(x,y) = w A (x) + w B (x,y) 


( 3 ) 


The first part, denoted by the subscript A, is an axisymmetric prebuckling displacement. 
The second part, denoted by B, is the infinitesimal nonaxi symmetric displacement that 
occurs at buckling. 


Prebuckling Equations 

The equilibrium equations which govern the axisymmetric prebuckling are obtained 
from substitution of the axisymmetric displacement in equations (3) into equations (1) and 
are found to be 


N = 0 


N xy A " 0 

- M x a + ^- N Xa W A + p = 0 


(4) 


where the primes denote total differentiation with respect to x. 

The appropriate set of boundary conditions is found similarly from equations (2) to 


be 




N x +N =0 
x A x 

or 

U A 

o 

ii 

or 

V A 



A 

A 

M + Ne = 0 
x A x 

or 

w A 

M ' + N w A = ° 
X A X A A 

or 

W A 


(5) 
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In equations (4) and (5), 



The first of equations (4) and (5) require that N x ^ = Constant = -N x . The second 
of equations (4) and (5) require that = Constant = 0 (no applied shear). From 

equations (6a) and (6b) together with N XA = -N x , Ny A can be written as 


N, ; 


3 

rM x E s A s z, 

S " 

TT7 J- 

E y( ] 


E xl 

( 1 " 

h 

< 

i 

1 

E X P 

1^1 - M X M yj 


N 

^ 6(x-jZ) 


E r A r 


i=i 


R 


W A " =T 




” M X M yj 


N-s 


(7) 


where 





E s A s 

d 


The equation determining the prebuckling deflection is now obtained from equa- 
tions (4), (6), and (7) as 


where 


01 W A + /3 2 w a + 


N 


03 + ^ 6(x-jZ)04 

j=l 


W A " ^5 = 0 


D x , E s I s , E s A s^s E x 
01 = + — — + — 


MxMy 


d E- 


“ M x Myj< 


( 8 ) 


(9a) 
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(9b) 


g2 = 5 2E ^ E ffe + N x 


e x R(i - P x £y) d 


03 =H 


i 


E x r2 (1 - Mx?y) 


e q a q 
* x + d 


(9c) 


E r A r 

04 =-¥ 

R z 


(9d) 


E-yZ-^X ** 

/35 = E X R (! “ MxMy)^ ‘ P 


(9e) 


- 2 


and I g = I 0 - A s z s . The boundary conditions are 


M Xa + N x 5 = 0 


or w A - 0 




m x A " N x w A = 0 or W A = 0 


( 10 ) 


The reciprocal relations M x Ey = j3 y E x and p x D y = have been employed to 

simplify the equations. 


Buckling Equations 

The equilibrium equations and boundary conditions which govern the buckling of a 
stiffened cylinder with discrete rings are obtained by substituting equations (3) into equa- 
tions (1) and (2). If only linear terms in the buckling displacements (subscript B 
terms) are retained, and if identities (4) are subtracted out, the following buckling equa- 
tions are obtained: 


Nx B,x + Nx ^B,y " ° 

%B,y + N *y B ,x = ° 

Ny B 

" Mx B,xx + Mx yB,xy " M yB,yy + “TT ' N x A w B,xx * N x B w A 
- N y A w B,yy “ 2N xy A w B,xy = 0 


(ID 
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with the boundary conditions 


and 


where 




N Xb = 0 

W 

ii 

o 

N xy B = 0 or 

V B = 0 

M X g = 0 or 

w B,x = 0 

x B,x 

+ n x a w B,x + n x b w A + N xy A w B,y = 0 

or 


Wg = 0 

J 


( 12 ) 


N. 


E, 


x 


B 1 - PxMyL 

E C A 


u B,x + w A w B,x + My(v B 


+ M 

>y r ) 


“S A S / . - \ 

+ ( u B,x + w A w B,x - z s w B,xx j 

% = 1 - p x M y t B ’ y + ^ + Mx ( UB ’ x + W A w B,x)j 

N / w 

+ ^ 6(x-j Z)E r A r fv B>y + ~ - ZrWB,yy) 

j=l 

N xy B = G xy( u B,y + v B>x + w A w B>y ) 


(13a) 


(13b) 


(13c) 


M, 


B 


r^Tr-(wB,xx + MyW B , yy ) + ^w B>xx 


E s A 0 z 


MxMy 

^( u B,x + w A w B,x - z s w B,xx) 


(13d) 


12 


M 


yB 


'-'p^Kyy + ^B.xx) 

N r 

+ Y 6(x-ji) E r I r w B>yy - E r A r z r (v Bjy + ^ - z r w B>yy ) 
. 1=1 L V - 


(13e) 


M. 


xy B 


and 


C S J S v 

2Dxy + “ + 2 6 ( x -^) & 

j=l 


I r - I Qr - Aj.Zj.2 


r J r 


w B,xy 


(13f) 


To solve the buckling equations (11) and apply boundary conditions (12), the prebuckling 
quantities (subscript A) must first be determined by solution of the prebuckling equa- 
tion (8) and boundary conditions (10). The solution of the prebuckling and buckling equa- 
tions is discussed in the next section. 


SOLUTIONS 

In this section solutions of the prebuckling and buckling equations are obtained for 
simply supported and clamped shells by use of the finite-difference method. The pre- 
buckling and buckling equations are solved by employing a technique in which the system 
of governing differential equations is changed to a system of second- order differential 
equations and then these equations are written in terms of finite differences. The 
resulting difference equations are then written in matrix form and solved by a modified 
Gaussian elimination technique. Details of the solutions follow. 

Solution of Prebuckling Equations 

The governing equation for the prebuckling displacement, equation (8), is a fourth- 
order, variable- coefficient, ordinary differential equation and is changed to a system of 
second- order differential equations by letting 


f M 

? = w A 


(14) 
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Now equation (8) can be written as 


01 % + 0 2 % + 


N 

0 3 + ^ 6(x-jZ)j3 4 

i=i 


W A 



(15) 


The terms /3j to /3g are defined by equations (9). Equations (14) and (15) are now the 
governing equations for the prebuckling displacements. 

When the derivatives occurring in equations (14) and (15) are replaced by finite 
differences, a set of simultaneous linear algebraic equations for the values of w^ and 
£ at discrete points along the length of the cylinder is obtained. The stations along the 
length of the cylinder are taken to be equally spaced and are numbered from i = 0 at one 
end of the cylinder (at x = 0) to i = k at the other end of the cylinder (at x = a). The 
following central-difference approximations for the first and second derivatives are used: 


(f'). _ f i+l ~ f i-l 
' h 2A 


ap"\ _ ^i+1 “ ^i + ^i-1 

V 1 li - n 


(16) 


where i indicates the station and A is the distance between adjacent stations. As a 
consequence of employing central differences, there are also stations corresponding to 
i = -1 and i = k + 1. When equations (16) are substituted, equations (14) and (15) for 
i = 1, 2, . . ., k - 1 become 


where 


k + i ~ 2 k + k- 


6 m&n 


h ■ ■— - A r — + Wi + 1^3 » ~ = 0 


W A i+ i ” 2 w Ai + w Aj_ ! 


-?i = ° 


(17) 


'iH = 1 

(i = 

i iH = o 

(i =£ 



( 18 ) 
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The subscript H denotes the difference station at which a ring is located and can be 
expressed as 


H = mj 


(i = 1, 2, . . N) (19) 


where m is the number of difference spaces A .between adjacent rings and N is the 
total number of rings on the cylinder as previously defined. Note in equations (17) that 
the finite-difference procedure requires that the ring properties be averaged (or smeared 
out) over one difference spacing A. At each of the stations, ^ and w A ^ are unknowns; 
hence, the total number of unknowns is 2(k + 1). Equations (17) represent 2(k - 1) 
equations. The boundary conditions yield the four additional equations that are needed to 
give 2(k + 1) equations and 2(k + 1) unknowns. 

From equations (10), with the assumption that the shell has the same boundary con- 
ditions at each end and is supported at the ends, the prebuckling boundary conditions at 
x = 0 and x = a are 


w A = 0 

m x a + N x 5 = 0 


or 


w A = 0 


( 20 ) 


By using equations (6), (14), and (16), these boundary conditions can be written in terms 
of finite differences at stations i = 0 and i = k (at the ends of the cylinder, x = 0 and 
x = a, respectively) as follows: 


WA 0 = ° 

WA k = 0 


( 21 ) 


For simply supported shells, 


-N 



(22a) 
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and for clamped shells, 


“ 2w Aj + 2w Aq = 0 

> 

A 2 4 k: 2 w Ak _ 1 + 2w Ak = 0 


In equations (22a), 


^6 



EsAgZ s 


(22b) 


(23) 


and is as defined by equation (9a). 

The governing difference equations (17) and boundary conditions (21) and (22) will 
now be solved by a matrix technique. Equations (17) can be written in matrix form as 


Z Ai l + LiZ A . + Z A . +1 = S (i = 1, 2, . . ., k- 1) (24) 

and equations (21) and (22) as 

KZ Al + LZ Aq = S (25) 

KZ A k -l + CZ A k = § W 

where 



(27a) 


„ *\\ (*■% ^4 6 1h) 

~er) 

-A 2 -2 


(27b) 
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s = < <3 


i2 % 


For simply supported boundary conditions, 


0 0 


0 0 


1 0 


0 1 


and for clamped boundary conditions, 


0 -2 


0 0 


A* 2 


0 1 


The set of matrix equations (24), (25), and (26) is solved by a Gaussian elimination 
technique which has been formalized in reference 15. This method proceeds as follows: 
The first of equations (24) together with equation (25) is solved for Zaj in terms of 
Za 2 ; this result is substituted into the next of equations (24) and Za 2 is found in terms 
of ZA 3 , and so on. Finally, the last of equations (24) together with equation (26) deter- 
mines Za^. j, and then all of the Za^ are calculated in reverse order. In the proce- 
dure just described, it follows from the first of equations (24) and equation (25) that 

Z Al = [Ll - L" 1 ^ 1 [za 2 + S - SL- 1 ] (29) 

Now the general result for Za^^ can be written in terms of Za^ + j as 


z Aj _ p i z A i+ i + Qi 
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From equation (29) and equation (30) for i = 1, Pi and Qj can be determined as 


4* 



( 31 ) 


The substitution of 


Z A . = Pi_iZ Ai + Qi_i into the general equation (24) gives 



(i = 2, 3, . . k - 1) (32) 


Thus 





(i = 2, 3, 


k - 1) 


(33) 


where Pi and Qj are given by equations (31). The recurrence- relation equations (33) 
provide all the Pj and up to Pjj-i and Qjj-1- By using equations (26) and (30), 
ZAfc_ i can b e found explicitly as 


where 



(34) 


Now Z Ak _ 2 > z Ak_ 3 ’ * • •> z Ai can be determined from equation (30). For this solution 
the only matrix inversions for the Z A ^ are of 2 x 2 matrices; therefore the solution is 
well adapted for rapid computer calculation. Once the Z Ai have been calculated, the 
prebuckling displacements w A at any station along the cylinder are known. 

Solution of Buckling Equations 

The buckling load of a stiffened orthotropic cylinder with discrete rings can be 
obtained by solving the previously developed buckling equations (11) along with boundary 
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conditions (12). In equations (11) and (12), the conditions of continuity around the cylin- 
der are satisfied if 


ug = U(x)cos 2Z 
R 

v B = V(x)sin ^ 

wg = W(x)cos 2Z 
R 

M Xb = M(x)cos ^ 


(35) 


where n, the number of circumferential waves, is an integer. The buckling equa- 
tions (11) may now be converted to ordinary differential equations with variable coeffi- 
cients. The variable M is included as an explicit dependent variable because the 
eighth- order system of differential equations will be expressed for convenience, as four 
second- order equations in u, v, w, and M. By making this change the finite- 
difference method of solution previously described for the prebuckling equation can be 
employed. 

In buckling equations (11), one additional step has to be performed in order to pre- 
vent the appearance of derivatives of x of higher order than two. A higher order 
derivative in wg will appear when the N X g term is differentiated in the first of equa- 
tions (11). This higher order derivative is eliminated as follows: Equation (13d) is 

solved for w B ^ and this result is substituted into equation (13a) to give N X g in 
terms of M X g and derivatives of x of less than second order. In performing this 
step and in substituting equations (35) and equation (7) for Ny^, the following variable- 
coefficient second- order ordinary differential equations for the buckling load are 
obtained: 


AiU" - A 2 U + A 3 V' + AjwXw" + (ajwa + A 4 )w' - A 2 w^W + A 5 m' = 0 (36a) 



N 


r n i 

v — ' 

a 3 u’ + a 6 v m - 

a 7 + ^ S(x-j£)Ag 
3=1 

V - A 3 vj' a w' - 

a 9 + a 10 w A + 6(x-jl)A n 

j=l 


(36b) 
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V 


(a 12 - a 13 wa)u' + 


r „ n 

Ag - A 14 wa + 2 j 6 (x-j^)A 11 
j=l 



N 


N 

+ 

- A l 5 + N x + a i 6 Wa - ^ 6 (x-jZ)A 1? 

w" + 

a 12 w a - a 13 w A w A _ 2, 6 ( x “^)' A 17 


3=1 


3=1 


+ 


N 

a 18 " a 19N x + a 20 w A + a 21 w A + ^ 6 (x-jZ)(A 2 2 + a 23 w a) 

3=1 


w - m" = 0 


(36c) 


-A 16 U' + D X W" - a 16WaW - a 2 4 W + M = 0 


(36d) 


where 


_ _ D x Eglg EgAgZg^ 

Uy — T T ' ’ ' 

1 - ju x ]Uy d d 

The coefficients Aj to A 24 are given in appendix D. Equation (36d) is obtained from 
equation (13d). 

Writing equations (36) in finite differences at each station along the length of the 
cylinder, using the same numbering system as in the prebuckling solution and difference 
equations (16), and then writing the resulting difference equations in matrix form results 
in the following matrix equation: 


where 


A iZBj_ 4 + BiZg i + CiZ Bi+1 - 0 


(i = 0 , 1, 2, . . k) (37) 


Z 


Bi 


r *\ 

Ui 

Vi 


Wi 


rj 


(38) 


and where the 4x4 coefficient matrices Aj, B^, and Cj are given in appendix D. In 
writing the matrix difference equations (37), the ring properties are again assumed to be 
averaged over one difference spacing A as in the prebuckling solution. Also, as in the 
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prebuckling solution, the boundary conditions give the additional equations needed in 
order to have the same number of unknowns as there are equations. 

From equations (12), with the assumption that the shell is supported at the ends, 
the boundary conditions on the buckling displacements are 


N Xb = ° 

or 

o 

u 

N xy B = 0 

or 

Vg = 0 


* 

w 

II 

o 


s 

w 

ii 

o 

or 

* 

w 

II 

o 


( 39 ) 


Equations (39) admit the possibility of eight different sets of boundary conditions (four 
simple support and four clamped). By substituting equations (35) and writing in finite- 
difference form, the eight sets of boundary conditions can be written in matrix form at 
stations i = 0 and i = k (at ends of cylinder, x = 0 and x = a, respectively) as 
follows: 


a O z B_i + B 0 Z B 0 + C 0 Z B! = 0 ( 4 °) 

A k z B k _i + ®k z B k + ^k z B k+ i = 0 

where Zg is as defined by equation (38), and A, B, and C are 4x4 coefficient 
matrices whose terms depend upon the selected boundary conditions and whose values 
are given in appendix D. 

In the 4x4 matrix coefficient expressions of equations (37), (40), and (41), the pre- 
buckling displacements and derivatives at each station are determined from the pre- 
buckling solution (eqs. (30) to (34), along with the appropriate boundary conditions). The 
matrix difference equations (37), (40), and (41) define a linear homogeneous system of 
equations which have a nontrivial solution only if the determinant of the coefficient matrix 

A 

equals zero. For a specified pressure the lowest value of N x for which this determi- 
nant is zero represents the buckling load of the cylinder. This value is found by solving 
the matrix difference equations by the modified Gaussian elimination method which was 
described in the prebuckling solution. That is, the first of equations (37) and equa- 
tion (40) are solved for Zg^ in terms of Zg^ to give 

z Bq = J]aoAo _1 Bo + B 0 ] [AoAo -1 Co - Cqjzgj (42) 
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The recurrence equation is of the form 


ZBj = p i z B i+1 (i = 0, 1, . . k) (43) 

and therefore Pq can be determined as 

P 0 = [AqAo^Bo + Bq) ^qAq^Cq - Co] - (44) 

Now substituting = Pi_iZg^ into the general equation (37) provides the result 

Pi = -[AiPi.i + Bj] "" 1 C i (i = 1, 2, . . ., k) (45) 

where Pq is given by equation (44). Finally, equations (41) and (43) provide the result 

pk p k-l p k + p kPk + ^kJ z B k+1 = 0 (46) 

Now if for a value of N x the determinant of the coefficient of Zg^^ in equation (46) 
is zero, then this zero is a zero of the coefficient matrix and this value of N x repre- 
sents a possible buckling load of the cylinder. The above solution requires only the 
inversion of 4 x 4 matrices for the determination of the buckling load and therefore is 
well suited for rapid calculation on a digital computer. 

The actual solution on the computer is done in the following manner: Initially 
an n, the number of circumferential waves, and an N x are selected (the initial values 
for n and N x were taken in most cases to be the values obtained from a classical 
solution, where w^ is constant). The prebuckling displacement and derivatives at each 
station are then calculated for the selected N x and the value of the determinant of equa- 
tion (46) is calculated; N x is then iterated until a zero of the determinant corresponding 
to the smallest value of N x is determined. Then n is changed and the above proce- 

A A 

dure is repeated to determine a new N x ; the lowest value of N x which is determined is 
then taken to be the buckling load. It should be noted that the presence of spurious singu- 
larities in the determinant of equation (46) sometimes make determination of the zero of 
the determinant difficult. This problem was resolved by employing the procedure 
described in reference 16. 

For all the computations a digital computer was used. Because the maximum 
order of the matrix was 4, a large number of finite-difference stations could be used 
without exceeding the storage capacity of the computer. For the calculations, the number 
of stations was increased until an increase in the number of stations did not affect the 
buckling load to within a specified accuracy. For the cylinders investigated, it was found 
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that approximately 200 to 300 stations along the length of the shell were required for 
proper convergence. 

In order to evaluate and compare the effect of discrete rings, the solutions of the 
prebuckling and buckling equations were obtained for the case in which the rings are 
assumed to be smeared out. This solution is presented in appendix B. It was also of 
interest to compare the buckling loads obtained herein with buckling loads obtained when 
classical prebuckling deformations (w^ = Constant) are assumed. For smeared rings 
the buckling solutions, with classical prebuckling deformations and simple- support 
boundary conditions assumed, are presented in reference 3. For discrete rings a 
Galerkin solution of the classical buckling equations is presented in appendix C. 


Buckling Shapes 

The shape of the buckling displacement along the length of the cylinder was deter- 
mined by employing the following procedure. The matrix of the coefficients of 
in equation (46) is defined in the following notation: 


^k p k-l p k + BP k + c k] 


a ll 

a 12 

a 21 

a 22 

a 31 

a 32 

a 41 

a 42 


a 13 

a 14 

a 23 

a 24 

a 33 

a 34 

a 43 

a 44 


(47) 


The solution for each of the four elements of the column vector the notation 

of equation (47), can be obtained as 



a 12 

a 13 

a 14 


u k+ i = 

a 22 

a 23 

a 24 

(48 a) 


a 42 

a 43 

a 44 



a ll 

a 13 

a 14 


v k+ i = - 

a 21 

a 23 

a 24 

(48b) 


a 41 

a 43 

a 44 
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(48c) 



an 

a 12 

a 14 

w k + l = 

a 21 

a 22 

a 24 


a 41 

a 42 

a 44 


a ll 

a 12 

a 13 

M k+1 = - 

a 21 

a 22 

a 23 


a 41 

a 42 

a 43 


(48d) 


Equations (48) determine the relative components of Zb^ + ^ and then the remaining 

are calculated in reverse order by using equation (43). Once the Zg^ have been 
calculated, the buckling displacements Wg at any station along the cylinder are known 
in terms of an arbitrary constant. Each buckling displacement is then divided by the 
absolute value of the maximum displacement to obtain the normalized buckling 

displacements 


w B 

Wg 

“max 


RESULTS AND DISCUSSION 

Because of the large number of geometric and physical parameters involved in this 
investigation, it is impractical to present results of a general nature. However, it is of 
value to present some computed results for large-diameter cylinders in order to study 
the influence of discrete ring stiffeners, prebuckling deformations, and eccentrically 
applied axial compressive load resultants (applied axial load and edge moment). Thus, 
computations of axial compressive buckling loads for boundary conditions of simple sup- 
port were made for two types of stiffened cylinders appropriate for large-diameter 
booster interstage structures: ring- and stringer- stiffened isotropic cylinders and ring- 

stiffened corrugated cylinders. In order to compare theory with the test results of ref- 
erence 7, computations of compressive buckling loads were also made for clamped inte- 
grally stiffened and Z- stiffened isotropic cylinders. 

The dimensions of the ring- and stringer- stiffened cylinders and of the ring- 
stiffened corrugated cylinders are given in figure 2 and the dimensions of the integral 
stiffened and Z- stiffened cylinders are given in reference 7. For the computations, the 
material in the cylindrical shell and the material in the stiffeners were taken to be 
identical, a value of 0.32 was assigned to Poisson's ratio for the stiffened isotropic 
cylinders, and the expressions given in reference 8 were used to compute the orthotropic 
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shell constants for the corrugated cylinder. Eccentricity effects of the stiffeners were 
studied by moving the stiffeners from the internal to the external surface of the shell. 

The results of these computations are presented in tables I to m and in figures 3 to 11. 
Details and discussion of the computed results follow. 

Ring- and Stringer- Stiffened and Corrugated Cylinders 

The ring- and stringer- stiff ened cylinders and the corrugated cylinders are sus- 
ceptible to buckling modes of panel instability or general instability, depending upon the 
number of rings on the cylinder and the eccentricity of the stiffeners. General instability 
is defined as the buckling mode in which the rings deform radially and the cylinder wall 
and rings buckle as a composite wall. Panel instability is defined as the buckling mode 
in which the rings have little or no radial buckling deformation and the cylinder buckles 
between rings. When discreteness of rings is accounted for, the mode of instability is 
determined by inspection of the buckling shape. 

For the calculations, simple- support boundary conditions of N X g = 0, vg = 0, 

Wg = 0, and M x ^ = 0 were assumed. The number of rings on the cylinders was varied, 
and nondimensional buckling loads N x /e x or N x /Et (where t is the effective thick- 

A s \ V 

ness of the isotropic cylinder wall, — p + tl were calculated by three separate theories: 

(a) Discrete rings, taking into account "exact" prebuckling deformations (eqs. (30) 

to (34) and eqs. (44) to (46)) 

(b) Smeared rings, taking into account "exact" prebuckling deformations (eq. (B2) 

and eqs. (44) to (46) with smeared- ring matrix coefficients of appendix D) 

(c) Discrete rings, assuming constant prebuckling deformations (eq. (C7)) 

The term "exact" prebuckling deformations refers to the deformations due to the loading 
and end- condition restraints that occur prior to buckling. Considering the prebuckling 
deformations to be a constant gives a buckling analysis which has a membrane state of 
stress prior to buckling (as done, for example, in refs. 1 to 5). This type of analysis 
will be referred to as "classical" buckling theory. Calculations were also made by 
employing the theory of reference 3 (eq. (15)), which is a classical buckling theory. For 
these calculations, general instability represents the case where the rings are smeared 
out, and panel instability is taken to be the buckling load which is obtained by considering 
a simply supported cylinder whose length is equal to the ring spacing. 

The nondimensional buckling loads for ring- and stringer- stiffened isotropic cylin- 
ders and for ring- stiffened corrugated cylinders are presented in tables I and II and are 
shown in figures 3, 4, and 5. The calculated values in the tables are given for the gov- 
erning mode of buckling, panel or general instability (panel- instability values are denoted 
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by a footnote). In the figures, the nondimensional buckling load N x /Et for ring- and 
stringer- stiffened cylinders and N x /E x for corrugated cylinders, is plotted against N, 
the number of rings on the cylinder. The buckling loads for N = 0 correspond to buck- 
ling of a cylinder with no rings on it. For the calculations presented in tables I and II 
and figures 3, 4, and 5, the stiffened cylinders are assumed to have the compressive load 
applied at the middle surface of the cylindrical shell (e = 0); therefore, the prebuckling 
state has a moment boundary condition equivalent to that of the buckling state 
( m x A = m x b = 0). 

The computed results for isotropic cylinders stiffened by internal rings and stringers 
are presented in figure 3. The curves shown in the figure are results obtained by using 
the classical theory of reference 3. The symbols denote results obtained from the pres- 
ent solution and the comparative solutions presented in appendices A and B. (A similar 
convention for curves and symbols is followed in figs. 4 and 5.) In figure 3, the coin- 
cidence of the curves and square symbols suggests that, in using classical theory, the 
differences resulting from treating the rings as smeared or discrete members are negli- 
gible. When the triangular symbols are compared with the solid curves, the effects of 
accounting for exact prebuckling deformations in general- instability solutions for 
smeared rings appear to be small. However, the large differences between the circular 
symbols and the curves indicate that when both discreteness and prebuckling deforma- 
tions are taken into account, as is done in the present theory, the predicted buckling loads 
are substantially lower than results obtained from classical solutions. Furthermore, 
although results of the classical solution indicate that a transition from panel to general 
instability will occur when the cylinder is stiffened by more than three rings, the present 
solution suggests that panel instability occurs in cylinders with as many as seven rings. 

A possible explanation of the large differences between classical buckling loads and 
mode shapes and those obtained from the present solution is associated with the assump- 
tion in the present calculations that the load is applied at the middle surface of the cylin- 
drical shell (e = 0). In the classical solution of reference 3, load is applied (implicitly) 
near the stringer- shell centroid. (See the discussion in the appendix of ref. 4.) Substan- 
tial effects of loading eccentricity (i.e., the effect of varying e) have been noted in refer- 
ence 6 for longitudinally stiffened cylinders, and the present results suggest that these 
effects persist in certain ring- and stringer- stiffened cylinders. 

For isotropic cylinders stiffened by external rings and stringers, figure 4 again 
indicates substantial differences between results obtained from classical solutions and 
results based on the present theory. As in the case of the internally stiffened cylinders, 
classical solutions with smeared rings and with discrete rings give virtually identical 
loads. The effect of accounting for prebuckling deformations in general- instability solu- 
tions for smeared rings is somewhat larger than that noted in figure 3 and is again 
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believed to reflect the loading eccentricity effects already mentioned. Note that the use 
of the present theory for cylinders with less than three rings yields panel-instability 
loads that are actually higher than classical buckling loads. Both theories predict tran- 
sition between panel instability and general instability at nearly the same ring spacing. 
The coincidence of the triangular and circular symbols on the figure is a consequence of 
the expected result that increasing the number of rings diminishes the effects of 
discreteness. 

For the ring- stiff ened corrugated cylinders studied, the computed results presented 
in table II and shown in figure 5 indicate that the effect of "exact” prebuckling deforma- 
tions on the buckling load is negligible. This result is expected because for an axially 
compressed cylinder, prebuckling deformations are a function of Poisson's ratio, and for 
a corrugated cylinder the value of the extensional circumferential Poisson's ratio Jly 
is very small. (See ref. 8.) The results for corrugated cylinders also show negligible 
effects due to discrete rings. 


Stringer- Stiffened Cylinders 

The influence of e, the eccentricity of the applied axial compressive load, on the 
buckling load is presented in table III and is shown in figure 6 for an isotropic cylinder 
stiffened by internal and external stringers (dimensions given in fig. 2(b)). For the cal- 
culations, simple- support boundary conditions N X g = 0, vg = 0, wg = 0, M x ^ = 0, 

W A = an d M Xa + N x e = 0 were assumed. In figure 6 the nondimensional buckling 
load N x /Et is plotted against the nondimensional distance e/h, where e is the radial 
distance from the cylinder's middle surface to a circumferential line on which the applied 
load N x acts and h is one-half the height of the stringer. (See fig. 2.) An eccen- 
trically applied compressive load is equivalent to an axial compressive load applied at 
the middle surface of the cylinder plus an axisymmetric end moment. The results pre- 
sented in table III and figure 6 show that the buckling loads for these cylinders are quite 
sensitive to the radial position at which the load is applied to the ends of the cylinder. 
Note that the maximum buckling loads of the configurations investigated can be as much 
as six times as great as loads obtained by considering alternate loading positions. This 
increase or decrease in load is an interesting phenomenon for it suggests that the design 
of end fittings may have an important influence on the strength of stiffened cylindrical 
structures such as launch-vehicle interstages. This effect, of course, does not occur 
with clamped cylinders. 

In figures 7 and 8, the experimental and theoretical buckling loads are shown for 
the integrally stiffened and Z.- stiffened test cylinders of reference 7. In these figures, 
the nondimensional compressive buckling load N x /Et has been plotted against the 
length- radius ratio a/R. Calculated curves are shown for the four sets of clamped 
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boundary conditions given by equations (39), for internally and externally stiffened cylin- 
ders. These curves were computed by employing the solutions given by equations (30) 
to (34) and equations (44) to (46), where, for the selected clamped boundary conditions, 
the proper terms of the selection matrix a are chosen (see appendix D). Also shown 
for comparison are the curves for the clamped solution of reference 7, which is a classi- 
cal buckling theory. The computed results in figures 7 and 8 show for the externally 
stiffened cylinder substantial changes in the buckling load due to taking into account 
"exact" prebuckling deformations. (Compare the classical-theory curves with the solid 
curves.) Substantial changes also occur as a result of considering different in-plane 
boundary conditions. (Compare the solid curve with the dashed curves.) For the inter- 
nally stiffened cylinder, there are only slight changes in the buckling load due to these 
effects. The agreement between theoretical and test results shown in figures 7 and 8 is 
generally good, and it is better when prebuckling deformations are considered than when 
they are neglected. 


Prebuckling and Buckling Shapes 

In figures 9 to 11, representative axial prebuckling and buckling shapes are shown 
for the simply supported stiffened isotropic cylinders previously presented in figures 3, 
4, and 6. In figures 9 to 11, the "exact" nondimensional prebuckling deformation w^/t 
or the normalized buckling deformation wg j | w B max | is plotted against the cylinder 
length x/ a. 

In figure 9, prebuckling and buckling shapes are shown for cylinders stiffened by 
external rings and stringers and by internal rings and stringers (all with four rings). 
Buckling results for these cylinders were presented in figures 3 and 4. The prebuckling 
shapes shown in figures 9(a) and 9(c) are quite different from the usual prebuckling shape 
of Foppl found in unstiffened cylinders (ref. 17). The buckling shapes illustrating 
general-instability and panel-instability modes of buckling, respectively, are shown in 
figures 9(b) and 9(d). 

Prebuckling and buckling shapes for the longitudinally stiffened cylinders of fig- 
ure 6 are shown in figures 10 and 11. The curves in figure 10 indicate that prebuckling 
shapes are sensitive to changes in the loading eccentricity parameter e/h and that 
large differences in prebuckling shapes exist between internally and externally stiffened 
cylinders with similar values of e/h. These results are reflected in the buckling loads 
presented in figure 6. 
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CONCLUDING REMARKS 


An analytical investigation has been presented of the buckling of eccentrically stiff- 
ened orthotropic cylinders, with discretely located ring stiffeners and "exact" prebuckling 
deformations included in the analysis. Nonlinear prebuckling equations and boundary 
conditions and linear buckling equations and boundary conditions have been derived by 
using energy principles, and solutions for simple- support and clamped boundary condi- 
tions have been obtained by the method of finite differences for cylinders loaded with any 
combination of axial compression and lateral pressure. Solutions were also obtained for 
"smeared" rings, and a Galerkin solution was obtained for the classical buckling problem 
(simple membrane prebuckling state) of a simply supported cylinder with discrete rings. 

Results of sample calculations have been shown for two types of axially com- 
pressed, simple supported, large- diameter, stiffened cylinders in order to illustrate the 
influence on the buckling load of prebuckling deformations, discrete rings, and eccen- 
tricity of applied axial compressive loads. In addition, the buckling loads of clamped 
longitudinally stiffened cylinders are presented in order to compare the present solutions 
with available test results. 

The calculations presented, which take into account "exact" prebuckling deforma- 
tions, indicate the following: 

1. The buckling load may be either substantially greater or less than that predicted 
by classical theory. 

2. The transition between panel- and general- instability modes of buckling is dif- 
ferent from that predicted by classical theory. 

3. The buckling load may be decreased for cylinders with few rings by considering 
the rings to be discretely located; thus, the common practice of "smearing" rings in 
design analysis can be unconservative. 

4. The buckling load may be substantially increased (as much as six times for the 
cylinders studied) or decreased, depending upon the radial position from the cylinder's 
middle surface at which the axial compressive load is applied at the ends of the cylinder; 
thus, the end-fitting design may have an important influence on the strength of stiffened- 
cylinder structures such as launch- vehicle interstages. 

5. However, the results for the ring- stiffened corrugated cylinders confirm the 
conclusion that effects due to considering prebuckling deformations and discrete rings in 
corrugated cylinders are negligible. 

6. The agreement between test results and theory is better than that previously 
achieved with theories that do not take into account prebuckling deformations. 
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The calculations show the sensitivity of stiffened cylinders to the investigated effects 
and, thus, demonstrate the need for calculations which accurately account for details of 
shell geometry and loading when investigating the stability of stiffened shell structures. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 13, 1967, 

124-08-06-26-23. 
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APPENDIX A 


DERIVATION OF GOVERNING EQUATIONS 

In this section the nonlinear equilibrium equations and boundary conditions for an 
eccentrically stiffened perfect orthotropic cylinder are derived by obtaining strain- 
energy expressions for the shell, stiffeners, and loadings and by applying the method of 
minimum potential energy. 


Strain-Displacement Relations 

For the coordinate system shown in figure 1, the Donnell type nonlinear strain- 
displacement relationships can be written as 


where 


e Xrp ~ e x ~ zw >xx 

e y T = e y “ zw, yy ? 

Ixy-p = I'xy ~ ^ zw ’xy 

> 

Yxy = u >y ^?x^>y 


e x ~ u >x + 2 W, x^ 


w 1 9 

e y = v, y + r + 2 w, y 


(Al) 


(A2) 


The quantities u, v, and w are the displacements of the middle surface of the ortho- 
tropic cylinder wall and a comma denotes partial differentiation with respect to the 
subscript. 


The stiffeners are assumed to behave as beam elements, 
strain-displacement relations can be written as 



and thus the stiffener 


(A3) 


where the subscripts r and s denote rings and stringers, respectively. Equa- 
tions (A3) specify that the displacements are such that the strain varies linearly across 
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the depth of the stiffener and satisfies compatibility of the displacements between the 
stiffener and the surface of the shell to which it is fastened. 


Energy of Stiffened Cylinder 

Orthotropic shell .- The strain energy of the orthotropic shell II C can be 
expressed in terms of resultant middle- surface forces and moments as 


p27rR r>a /_ _ 

Hc = 2 J 0 J 0 ( N x e x + N : 


xy^xy Ny e y - ^x^’xx ^^xy^xy “ ®^y^ ’ yy^ (A4) 


where (see ref. 14) 


N x = 

x ( c 



ii 

E y L 

1 " MxMy 

21 

$ 

11 

= G xy>xy 

M X = 

D x 

1 - y 


D v 

M v = 

y 

J 

1 “ Mx^y 

Mxy = 

- D w 
xy >xy 


e x + M v e 


r) 


(A 5) 


and e x , e-y, and y^y are given by equations (A2). 

Rings .- If the rings are considered to be equally spaced and discretely located 
along the length of the cylinder, the strain energy n r of N rings can be expressed as 


"r = | I i 0 2 ’ R (i’ Ar E W I > r + 


(A6) 


where dA r is an element of the cross-sectional area of the ring A r , G r J r is its 
twisting stiffness, and l is the ring spacing. In equation (A6) it is assumed the total 
circumferential arc length of the ring is approximately the same as that of the middle 
surface of the shell. 
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In the energy expression, the first term is the expression for the energy associ- 
ated with bending and extension of a beam and the second term is inserted as an approxi- 
mation to the energy due to twisting of the ring. The approximation for energy due to 
twisting results from assuming that the ring twists in such a fashion that its angle of 
twist is equal to the angle of twist of the shell. 

When the ring strain-displacement relation (eqs. (A3)) is substituted into equa- 
tion (A6) and the integrations are performed over the ring area, the strain energy for 
N rings becomes 

N „ „ _ 


where z r 
shell, I Qr 
5(x-jZ) is 


i V C*™ P d 7 2 — 2 

n r = g 2 _, J 0 Jq \®* r A r 6y “ 2E r A r z r w,yy€y + E r I 0i> w,yy 

+ G r J r w, xy 2 )8(x- j l )dx dy (A7) 

is the distance from the centroid of the ring to the middle surface of the 
is the moment of inertia of a ring about the middle surface of the shell, and 
a Dirac delta function defined so that 



f(x)6(x-jZ)dx = f(j Z) 


(x = jZ) 


6(x-jZ) = 0 


(x * j Z) 


(A8) 


S tringers .- If the stringers are assumed to be closely spaced so that their elastic 
properties can be averaged over the stringer spacing, the strain energy of the stringers 
n s can be taken as 



P 27rR r a /EsAs e 2 
Jo \ d £x 


2EgA s z s 

5 e x w ’xx + 


Eol, 


S A Oc 


' w ’xx + 



dx dy 


(A9) 


where d is the stringer spacing and the subscript s is used to denote stringer 
properties. 

Applied loads .- The potential energy II associated with applied loads of con- 
stant pressure p (positive for external pressure) and an externally applied load 

At 

resultant N x (positive in compression) is 

p2jrR pa p2irR ~ . _ a 

IIt = \ \ pw dx dy + \ N x (u - 5w ) dy (A10) 

J 0 J 0 J 0 y o 
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where e is the distance from the middle surface of the orthotropic shell to the line on 
which the load resultant N x acts. The distance e is introduced so that loadings of 
applied edge moments and axial loads may be represented by a statically equivalent force 
system. 

The total potential energy of the stiffened cylinder II is the sum of the energies 
given by equations (A4), (A7), (A9), and (A10): 

n - n^, + Hj» + rig + rijj (Ail) 

Equilibrium Equations and Boundary Conditions 

The nonlinear equilibrium equations and boundary conditions are obtained from 
equation (All) by applying the principle of minimum potential energy, which requires the 
vanishing of the first variation (611 = 0). After integration by parts and proper 
regrouping, the nonlinear equilibrium equations for an eccentrically stiffened orthotropic 
cylinder with discrete rings obtained in the above manner are 


N x,x + N xy,y - 0 


N y,y + N xy,x - 0 


> (A12) 


"^x,xx + Mxy,xy “ ^y,yy + “p" " ^x w ’xx ” ^y w >yy ” ^-^xy w >xy + P _ ® ^ 


with the following boundary conditions to be satisfied at each end of the cylinder: 


N x + N x =0 


Nxy'O 


M x + N x e = 0 


or u = 0 

or v = 0 


or w, x = 0 


Mx,x - Mjjy y + N x w, x + NjjyW^ - 0 or w 0 
The stress and moment resultants in equations (A12) and (A13) are 


N x = 


E, 


1 ~ /^x^ y 


u >x + 2 W,X + M y( v ’y + R + 2 W, y ) + d ( U ’ x + 2 W ’ X ~ ^s w, xx) (A 14 a) 


(A13) 
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N y - rnfev [> + 1 + l w, y 2 + + |'y.x 2 | 


N 


M x = - 


+ ^ 6 ( x ~i^ E r A r( v >y + ^ + 2 W >y “ z r w ?yy) 

3=1 

^xy = ®xy( u >y + v >x + w >x w >y) 

Dy EsIo s E s A s z s / j 2\| 

[i - MxM ~( w,xx + M y Wj yy) + — T" w,xx d v u,x + 2 WjX J 


(A14b) 


(A14c) 


(A14d) 


M, r = - 


D,; 


/i x jLt y ( w, yy + MxW,xx ) 


n r 

^ 6(x-jZ) E r I 0r w,yy - E r A r z r ^v,y + ^ + ' 2 w >y 

j=l L 


M xy “ 


2D 


N 


xy 


+ S ^ S + ^ 6(x-jZ)G r Jj 

j=l 


w, 


xy 


(A14e) 

(A14f) 


Equations (A12) and (A13) appear in the body of the paper as equations (1) and (2). 

f 
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SOLUTION OF PREBUCKLING AND BUCKLING EQUATIONS 
FOR SMEARED-OUT RINGS 

Prebuckling Solution 

It is of interest to present the solution of the prebuckling equation for the case in 
which the rings are assumed to be smeared out over the ring spacing in order to com- 
pare it with the discrete- ring case. When the rings are smeared out, the discrete- ring 
prebuckling governing equation ( 8 ) becomes the following constant-coefficient differential 
equation: 


+ ^2 W A + ^3 + "^) W A " % " 0 ( B1 ) 

where l is the ring spacing and to /3g are defined by equations (9). Equa- 
tion (Bl) has the following analytical solution that satisfies boundary conditions (10): 


= yi sin 0 !]X sinh a^x + 72 s ^ n a l x cos h a 2 x + 73 cos a l x cos h 0L 2 X 


where 


+ 7/4 cos oi jx sinh a 2 x + 


05 



a l = 



a 2 


'£ 02 _ 

2 40i 





J 


(B2) 


(B3) 


For simply supported shells, 


%( a 2 2 ~ a l 2 ) | % 


(B4a) 
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(cos ofja - cosh a 2 a)< 


V 2 


Linh Q! 0 a 

I 

06 [ 05 ( a 2 2 " "l 2 ) 

~\ 

2a'ja'2|3j5 sin a^a 

01 . 04 

. 04 { 

l 

ft 3 + T J 

03+ T 

J 


2«i« 2 (cos “i 3 - sinh a 2 a )^ + ( sin a l a cos ^ a 2 a )^ 


(cos aja - cosh a^a^ 


V 4 = 


\. 

Sin Q!^ 

06 05(“2 - “1 ) 

sinh a! 2 a 

ftl „ 04 

„ % [ 

1 

V 

L 3 T 

03 + T 

y 


r l = 


^“1^2 J( cos a i a sinh Q! 2 a ) 2 + ( sin a l a cosh a 2 a )^j 

For clamped shells, 

^5 ( a 2 2 + a l 2 ) sin a l a sin ^ a 2 a ” Q! i Q: 2( s ^ n ^^ Q! 2 a + s i n ^ a j a ) 


h 

h + T 


2.2 2 . . 2 

Of 2 sin Q!j a - oij sinh o^a 


V2 = 


j 3 g (cos a!ja - cosh a^a^a^sin a^a - &2 a \ sin ^ a 2 a ) 


0 4 

*3 + T 


2 2 2 2 
0*2 sin Q!j a - sinh a^a 


r 4 = 


(cos a!ja - cosh atga^o! j^sinh a^a - a \ a 2 s ^ n “I s ) 


04 

* + f 


2.2 2 . u 2 
ag sin a^a - aj sinh o^a 


(B4b) 


(B4c) 


(B5a) 


(B5b) 


(B5c) 


In equations (B4), /3g is defined by equation (23). Similar types of solutions for pre- 

buckling displacements are reported in references 6, 9, and 10. 

Buckling Solution 

When the rings are assumed to be smeared out over the ring spacing, the buckling 
equations (11) and boundary conditions (12) are still applicable except that the stress and 
moment resultants defined by equations (13) are defined as follows: 
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„ Ey T W B ~ / » 71 E r A r( W B - \ 

%B " l - MxMyL B,y + 'r“ + Mx ( Ub » x + w A w B,x)| + — — ( v B,y + ~ z r w B,yyj 


M- 


YB 


D El 

l - MxMy ( WB »yy + MxWB » xx ) + “T^ w B,yy 


Ep A 


Mr/ 

— l V B, 


Wt 


y + __ _ ZrWB 


»yy) 


> (B6) 


m xyb “ ( 2D x y + “V + “T^j^xy 


and N X g, M x ^, and N X y B are as previously defined by equations (13) and l is the 
ring spacing. Similar equations for stiffened cylinders are presented in reference 4 
except that an isotropic shell is used in reference 4 where an orthotropic shell is con- 
sidered herein. To solve the smeared-ring buckling equations (11) for the boundary con- 
ditions considered herein, the finite -difference method is used. The procedure is exactly 
the same as previously described for the discrete -ring case, and equations (35) to (46) 
are the same except for the following three changes: the prebuckling displacements and 

derivatives are given by the analytic solution of equation (B2); the summation and the 

N 


Dirac delta function 


£ Ofc-j*) 
3=1 


are replaced by 


l/l in the buckling equations (36) 


(Aj to Ag 4 are defined the same); and the matrix coefficients A}, Bj, and Cj of 
equations (37) are defined in appendix D. 


For the numerical computations, a computer program similar to the program for 
discrete rings is used with the smeared-ring matrix coefficients Aj, Bj, and Cj, and 
with the prebuckling displacements given by equation (B2). The same procedure of cal- 

A 

culating the buckling load N x is used. 
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CLASSICAL BUCKLING SOLUTION FOR DISCRETE RINGS 


When the classical prebuckling displacement state is assumed (that is, 

W A = Constant), the classical-theory buckling equations for a stiffened orthotropic cylin- 
der with discrete rings can be determined from equations (11) as 


N x + N™ = 0 
x B,x ^Bjy 


N YB,y + Nx YB,x " ° 


■'v 


N 


YB 


" Mx B,xx + Mx YB,xy " M YB,yy + ~R~ + N x w B,xx + P Rw B,yy = 0 
with boundary conditions of 


o 

II 

or 

o 

ii 

PQ 

3 

Nxy B «0 

or 

- v B = o 

m *b = 0 

or 

w B,x = 

Mx B,x ' M *y B> y ' ” xW B,x = 0 

or 

w B = 0 


J 


where 


N< 


X B 


Mx^y 


/ W B \ E S A S _ 

u B,x + Py( v B,y + j + J~( u B,x “ z s w B,xx) 


Ny B = i ( v B,y + ^ + Px u B,x) + J 6 ( x -30E r A r (v B> 

x Y j=l 

N xy B = G xy( u B,y + v B,x) 


W B - 

y + R - Z r W Bj yy 


(Cl) 


(C2) 


(C3a) 


(C3b) 


(C3c) 
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Du / v E s ! s Eg A z 

( w B,xx + My w B,yy) + — wb,xx ( u B,x " z s w B,xx) 


1 - Mx^y 


{ IX 

7 ~TT ( w B,yy + M x w B>xx ) + Z 6(x - ]l) E r I r w B,j 
A " M X My j=l L 

- E r A r z r ^VB > y + - z r w B,yyj\ 


Mxyg = ^^xy + , + ^ $( x- j OG r J r Wg^xy (C3I 

L i=i _ 

and N x is the compressive load, Ny^ = -pR, and Nxy^ = 0. 

To obtain a solution for the buckling load which satisfies boundary conditions (C2) 
of classical simple support ^vb = N X g = M X g = wb = oj, the buckling displacements are 
taken as 

~\ 

OO 

u B - cos f l c q cos 3S 


Vg = sin 


OO 

Z b q sil 


Wg = COS 


Z a q siI 


where n is the number of circumferential waves and the coefficients aq, bq, and Cq 
are to be determined. 

Use of the Galerkin method and a procedure similar to that of reference 11 yields 
the following set of equations for the unknown coefficients aq, bq, and Cq: 
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Q x (u,v,w)cos ^ cos dx dy = 0 

Q v (u,v,w)sin ^ sin LEL dx dy = 0 
•* k a 

Q z (u,v,w)cos ^ sin dx dy = 0 


*\ 




J 


/s = 0, 1, 2, . . . 
\r = 1, 2, 3, . . . 


(C5) 


where Q x , Qy, and Q z are the result of substituting the series for the buckling dis- 
placements u, v, and w, equations (C4), into the governing equations (Cl). 

From equations (C5) the resulting system of equations is determined as 


■Fjc r + F 2 b r + Fga r = 0 


2Fc 

F 2 c r - F 4 b r - ~ 


N 


R q=l 


N 


- Fe^r-T 1 I sin iT?I I 

j=l q=l 


N 


2F 7 

R 

x2 ,„.2 


F 3 c r “ F 6 b r - V- 2 Sin iT?T Z b( l Sin 5??T 

q=i 


“ p - N x (?) - P*(£) 

_ 2 ( F 9 + F _ lo) y sin^L y a^ sin _aiQ_ 

a L N+R ^ N+l 


2Fn 


a 


N oo 

V ^cos-llL y * SUL cos ^1^=0 

L a N+l A 4 a N + l 


3=1 


q=l 


> (r = 1, 2, . . .) (C6) 


where Fl to Fjj are defined in appendix D. Equations (C6) give the criteria for 
buckling if the determinant of the coefficients vanishes, that is. 


APPENDIX C 


By combining equations (C6), the following determinant of the coefficients a^. can 
be obtained: 


*rq 


Fg + (F 9 + F 10 )$ r + F n t, 


I — 

+ F F 2 (F 6 + F 7 $ r ) - Fg(F 4 + F 5 $ r ) 
3 Fj(P 4 + F 5 4 r ) - F 2 2 

+ (Fg + F 7 * r ) F 3 F 2 - F l( F 6 + F 7*r) 
^ F l( F 4 + r 5^rJ - f 2 2 


- »x(¥) - P R (g) 


= 0 


6 rq + F 9$rq + F ll*rq 

(r = 1, 2, . . q = 1, 2, . . .) (C7) 


where 


* I I sln2 

1-1 


rjrj 

N + 1 





rq 


N 

=11 

i=i 


rn qn 
a a 


cos 


N + 1 


cos 


_aiL 

N + 1 




J 


(C8) 


Examination of the determinant given by equation (C7) reveals that the terms with 
r and q odd are not coupled with the terms with r and q even. Thus, the deter- 
minant of equation (C7) may be separated into two determinants, one for r and q odd 
and one for r and q even. The determinant for r and q odd represents a sym- 
metrical w-deformation and the determinant for r and q even represents an antisym- 
metrical w-deformation. The form of these two determinants is identical to the deter- 
minant of equation (C7) except for the mentioned restrictions on r and q. Each 
determinant given by equation (C7) (r and q odd,, or r and q even) is in the form of 
a standard eigenvalue problem and is solved on a digital computer. When computing the 
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buckling load, each determinant is minimized numerically for integral values Qf n, the 
number of circumferential waves in the buckle pattern, and the buckling load is taken to 
be the absolute minimum obtained from the r- and- q- odd determinant and the r- and- q- even 
determinant. The size of the determinants (r) is determined by requiring that the buck- 
ling load converge to a desired accuracy. For the calculations presented herein, the size 
of the determinants was increased until the buckling load between a determinant of a 
specified order and a determinant of one less order did not vary to within a specified 
accuracy. For the cylinders considered herein, a maximum of 10 terms was required 
for convergence of the series given by equations (C4). 
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FORMULAS FOR COEFFICIENTS 

The coefficients Aj, A£, . . A24 appearing in equations (36) are as follows 

E y 


Ai = 


E S A S ( D x 




W 


1 “ Mx^y D x d yl “ l^x^y ^ 

2 


a 2 - G xy(p) 

Ag. f .- W 

\1 ~ MxMy 7 R 

^x^y D x ^yE s AgZ s /n\^ 


A 4 = 


" Mx^y) D x (l " Mx^y)^ 




E s A s z s 
A 5 = S" 


D x d 


A 6 _ G xy 


A 7 = 


Ey -.2 


1 - 


MxMy 


(I) 


A g - E r A r 


m 


A 

A » = ^27 — 


R (l Mx^y) 


A 10 " G xy R 


A ll = 


E r A r n 

1 

R 

R 
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EyP X 


R(1 - Mx^y) 


E x + E s A s 
1 - Mx^y d 


E xMy _n 
- /ixMy E 


M ^ +2D ^)® 2 


E s A s z s 
A 16 = — d 


a 17 “ G r J r(R) 


Aio = -z * 

J-o 1 - a ii 


/Jlf + 

W R 2 (l - P x P y ) 


a _ ^y^x < 
19 ” E x(! - M x My)' 


Aon = — 


E yMx E s A sZs /n \ 2 _ E yffx 
E x (l - Mx^y)d ^ " ^x^y) 


21 E X R(1 - M x My) ^ R > 


/ o\ / n \4 2E r A r z r n 2 EyAj» 

A 22 = (e* + E r A r $ r 2 )(|) + -J + — r 


A 23 “ 


E r A r n 2 


Ao/i = 


®x^v / n> 


24 l . ^/iy' 
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The coefficient matrices Aj, Bj, and Cj appearing in equation ( 37 ) are as 
follows. 

For discrete rings: 


Al 

aa 3 

2 

1 / » A m \ . ^ 

l w A - 2 Wa ) Ai " 2 

1 

AA 5 

2 

AA 3 

2 

a 6 

AA 3 w A 

2 


0 

f-( A 13 w A " A 12 ) 

0 

_A 15 + N x + A 16 W A + |(A 13 w; 

A 17 |£fi+ 11 H_ 6 _ 

A L 4 lH 4 J 

k ~ A 12 ) w A 

-1 

AAi 6 

2 

0 

aa 16 , 

D x + 2 W A 


0 


-2Ai - A 2 A 3 


"( 2A 1 + A 2 A 2 )w A 


-2Ag - A Arj - AA g 6 iH 


-A A 9 - A A 1Q w A - AA 11 8 iH 


B i = 


A 2 A 9 - A 2 A 14 w a 
+ AA ll 6 iH 


2A 4 5 + A 2 A 4 g - (2 + A 2 A j 9 )n x 
+ (a 2 A 2 o - 2A 16 )w A + A 2 A 21 w a 
/2Ai7 


\ A 


,+ AA 22 + AA 23 w A 6 iH 


-2D X - A A 24 
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A i 

aa 3 

2 

/ t A T T \ A ^ ^4 

( W A + 2‘ w A.j A l + 2 

aa 5 

2 

AA3 

2 

A 6 

aa 3 w A 

2 

0 

J ( A 12 - a 13 w A.) 

0 

- a 1 5 +N X + A 16 wa 

+ f( A 12 " A 13 w a) w A 
+ *17 . 6ih _ 

-1 

aa 16 

2 

0 

_ AAi 6 WA 
x 2 

0 


where the subscript H of the Kronecker delta is as defined by equation (19). 
For smeared rings: 


— 




A 1 

A A 3 
2 

( W A - f w A.) a 1 " f a 4 

aa 5 

2 

AA 3 

2 

A 6 

AA 3 , 

,^ w A 

0 

f ( a 13 w A “ a 12) 

0 

-A 15 + N x + A 16 wa 
+ 2 ( A 13 w A - a 12 )wa - * 7 

-1 

AA 16 

2 

0 

_ AA 16 wa 

x + 2 

0 
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-2A^ - A^A2 

0 

~[ 2A 1 + A 2 A 2 )w^ 

0 

0 

-2A 6 - A 2 A 7 - 

-a 2 a 9 - a 2 a 10 wa - t n 

0 

0 

A 2 A 9 - A 2 A 14 wa 
A^A n 
+ l 

2A 15 + A 2 A 18 - (2 + A 2 A 19 )n x 
+ (a 2 A 20 - 2A 16 )wa 
+ A 2 %wa+2 ^ +a 2M 

a2a 23 

+ 1 W A 

2 

0 

0 

-2D X - A 2 A 24 

A 2 







A 1 

aa 3 

2 

/ T A T T \ A ^ A4 

(*A+2 w A)Al+ 2 

aa 5 

2 

aa 3 

2 

A 6 

aa 3 w A 

2 

0 

f ( A 12 - a 13 w A.) 

0 

_a 15 + N x + Ai 6 wa 
- f ( a 13WA - a 12)wX - 

-1 

aa 16 

2 

0 

_ A A i6 wA 
x 2 

0 
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The coefficient matrices Aq^, Bq^, and Cq^) appearing in equations (40) 
and (41) are as follows: 


'O(k) 


= a 


aa 13 


0 

AAxe 


A6 

2 

0 

0 


aa 13 * A 
2 W A - A 16 


0 

0 

0 


0 

0 

0 

0 


B 0(k) " a 


0 

a 2 a 14 

a2a 12 + A 16 

0 


1 - an 

0 

0 

0 

" AA 10 

0 

-AA 10 w A 

0 

+ 

0 

1 " “22 

0 

0 

0 

0 

0 

0 


0 

0 

1 

0 

0 

0 

-2D X - A 2 A 24 

A 2 


0 

0 

0 

1 

Ql 

1 


where 


'O(k) 


= a 


AA 


13 


2 

0 

0 

aa 16 


a = 


“11 

0 

0 

0 


A6 

2 

0 

0 


0 

“22 

0 

0 


- A le 


0 

0 

2Dx 


0 

0 

0 

0 


0 

0 

0 

“44 


0 

0 

0 

0 
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The elements <244, and «44 of the selection matrix a take on the value 1 or 0 

depending upon the prescribed boundary conditions. The selection matrix a is used to 
select the prescribed boundary conditions in the following manner: if «22> or 

(244 equals 1, then N X g = 0, Njjyg = 0, or wj^ x = 0, respectively; and if an, <*22 , 
or 544 equals 0, then Ug =0, Vg = 0, or M X g = 0, respectively. For the simple- 
support boundary conditions considered herein (N X g =0, Vg = 0, M Xb = 0), 5 44 = 1, 

<*22 = 0> and <244 = 0, respectively. 


F l = 


coefficients 

F l> 

f 2 » • • •> F 11 

( E x + 

E S A 

sVnrf + G / 

\1 - MxMy 

d 

yUJ +<Jx y( 

/ E x p y ^ 

G ) 

rw n 

yl - Mx^y 

°xyl 

a R 

E x^y 

£1, 

E s A sZs/r7rf 

R(1 " MxMy) 

-f 

a 

d l a,/ 


F 4 = 


1 - 


+Gx y (?) 


f 5 _ E r A rl 


z 

(I) 


Fe = 


E y n 


R^(l " Mx^y) 


F7 = E r A r 




F = ( ° x + EsIs + E s A s^s 2 V r7r\ 4 + / 2D x M y + 2D + ^s^s\ /ii\ 2 /rffY 

^ ^ ^ /\ a / yl - Mx/Xy d J \R/ \ a / 

+ Bz /nf 

1 - M X P y \R/ 


R 2 (l - M-x/^y) 
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TABLE I.- BUCKLING CALCULATIONS FOR SIMPLY SUPPORTED 
AXIALLY COMPRESSED RING- AND STRINGER- STIFFENED 
ISOTROPIC CYLINDERS 


Number of 
rings 

Discrete rings, 
"exact" 
prebuckling 
deformations 

Smeared rings, 
"exact" 
prebuckling 
deformations 

Discrete rings, 
"classical" 
theory 

"Classical” 
theory of 
reference 3 

N 

N x /Et 

n 

N x /Et 

n 

Nx/Et 

n 

Nx/Et 

n 


Rings external; stringers external 




0 

0.00085 

16 

0.00085 

16 

0.00058 

19 

0.00058 

19 

1 

a . 00238 

18 

.00314 

6 

a . 00134 

25 

a . 00134 

25 

2 

a .00284 

19 

.00323 

7 

a . 00241 

26 

a . 00241 

26 

3 

.00320 

7 

.00351 

7 

a . 00361 

18 

a . 00361 

18 

4 

.00353 

6 

.00368 

6 

.00447 

7 

.00447 

7 

5 

.00376 

6 

.00382 

6 

.00463 

7 

.00463 

7 

6 

.00391 

6 

.00395 

6 

.00478 

7 

.00478 

7 

7 

.00403 

6 

.00415 

5 

.00493 

5 

.00493 

5 



Rings internal; stringers internal 




0 

0.00017 

16 

0.:00017 

16 

0.00022 

15 

0.00022 

15 

1 

a .00041 

27 

.00211 

9 

a . 00064 

17 

a . 00064 

17 

2 

a .00068 

62 

.00242 

8 

a . 00139 

18 

a. 00139 

18 

3 - 

a .00095 

89 

.00262 

8 

a . 00244 

19 

a.00244 

19 

4 

a . 00129 

112 

.00281 

8 

.00299 

7 

.00299 

7 

5 

a . 00162 

138 

.00293 

8 

m .00309 

7 

.00309 

7 

6 

a . 00214 

162 

.00297 

8 

.00320 

7 

.00320 

7 

7 

a . 00263 

183 

.00304 

8 

.00329 

8 

.00329 

8 


a Buckling in panel- instability mode. 


55 



TABLE H.- BUCKLING CALCULATIONS FOR SIMPLY SUPPORTED AXIALLY 
COMPRESSED RING- STIFFENED CORRUGATED CYLINDERS 


Number of 
rings 

Discrete rings, 
"exact" 
prebuckling 
deformations 

Smeared rings, 
"exact" 
prebuckling 
deformations 

Discrete rings, 
"classical" 
theory 

"Classical" 
theory of 
reference 3 

N 

N x /e x 

n 

n x /e x 

n 

1 n x /e x 

n 

1 n x /e x 

n 




Rings internal 





0 

0.00050 

0 

0.00050 

0 

0.00050 

0 

0.00050 

0 

1 

.00096 

6 

.00098 

6 

.00097 

6 

.00098 

6 

2 

.00115 

6 

.00116 

6 

.00116 

6 

.00116 

6 

3 

.00133 

6 

.00133 

6 

.00133 

6 

.00133 

6 

4 

.00149 

6 

.00150 

6 

.00150 

6 

.00150 

6 

5 

.00166 

6 

.00166 

6 

.00166 

6 

.00166 

6 

6 

.00182 

6 

.00182 

6 

.00182 

6 

.00182 

6 

7 

.00198 

6 

.00198 

6 

.00198 

6 

.00198 

6 




Rings external 





0 

0.00050 

0 

0.00050 

0 

0.00050 

0 

0.00050 

0 

1 

a .00197 

0 

.00290 

0 

a .00197 

0 

a .00198 

0 

2 

.00326 

0 

.00337 

0 

.00336 

0 

.00337 

0 

3 

.00380 

0 

.00383 

0 

.00383 

0 

.00383 

0 

4 

.00428 

0 

.00429 

0 

.00429 

0 

.00429 

0 

5 

.00474 

0 

.00475 

0 

.00475 

0 

.00475 

0 

6 

.00520 

0 

.00521 

0 

.00522 

0 

.00522 

0 

7 

.00567 

9 

0 

.00567 

0 

.00568 

0 

.00568 

0 


a Buckling in panel- instability mode. 



TABLE IE.- INFLUENCE OF ECCENTRICALLY APPLIED COMPRESSIVE 
LOADS ON BUCKLING LOADS OF SIMPLY SUPPORTED 
STRINGER- STIFFENED ISOTROPIC CYLINDERS 


e/h 

N x /Et 

n 

Stringers internal 

1.0 

0.000111 

17 

0 

.000171 

16 

-1.0 

.000298 

14 

-2.0 

.000489 

13 

-3.0 

.000618 

13 

-4.0 

.000649 

16 

-5.0 

.000630 

16 

-6.0 

.000599 

17 

Stringers external 

2.0 

0.000220 

26 

1.5 

.000312 

24 

1.0 

.000494 

20 

.5 

.000740 

17 

0 

.000851 

16 

-.5 

.000841 

19 

-1.0 

.000724 

32 

-1.5 

.000628 

36 

-2.0 

.000557 

41 
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0.10 

(0.25) 


TT 


0.10 

(0.25) 


2.00 

(5.08) 


4.00 

(10.16) 


Ring section 


(a) Ring- stiffened corrugated cylinders, a = 200 in. (508 cm). 



i 

^ 2.50 _ 
(6.35) 

Stringer 


Li 3 

, 3 


0.10 

(0.25) 

1 

q I | Stringer section 


R - 200 (508) 


0.10 

(0.25) 


2h = 2.00 
, (5.08) 


L, 0.90 

1 (2.29) 


(b) Ring- and stringer- stiffened isotropic cylinder. 

a = 200 in. (508 cm). (See fig. 2(a) for ring section.) 

Figure 2.- Dimensions of stiffened cylinder. Dimensions are in inches. 
(Parenthetical dimensions in cm.) 
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.005 


N x 

If 


.004 k 


O Discrete rings ("exact" prebuckling deformations) 

A Smeared-out rings ("exact" prebuckling deformations) 
□ Discrete rings ("classical" theory) 

— Panel instability 

— General instability J " Classical " theor y of reference 3 
Solid points denote panel instability 

Open points denote general instability 



N 


Figure 3.- Buckling predictions for simply supported axially compressed isotropic 
cylinders stiffened by internal rings and internal stringers. 
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0 

Figure 4.- 




O Discrete rings ("exact" prebuckling deformations) 

A Smeared-out rings ("exact" prebuckling deformations) 
□ Discrete rings ("classical" theory) 

Panel instability 

— General instability J " Classical " theor V of reference 3 
Solid points denote panel instability 
Open points denote general instability 


j i i i i i i 

1 2 3 4 5 6 7 


N 

Buckling predictions for simply supported axially compressed isotropic 
cylinders stiffened by external rings and external stringers. 


.006 r 


.005 


.004 


.003 


.002 


.001 


O Discrete rings ("exact" prebuckling deformations) 

A Smeared-out rings ("exact" prebuckling deformations) 
□ Discrete rings ("classical" theory) 






Figure 5.- Buckling predictions for simply supported axially compressed ring- 

corrugated cylinders. 


stiffened 
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Figure 6.- Influence of eccentrically applied compressive loads on the buckling loads of simply supported stringer 

stiffened isotropic cylinders. 





(a) Prebuckling shape; external rings and external stringers. 



0 .5 1.0 


x 
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(b) General instability; external rings and external stringers. 

Figure 9.- Prebuckling and buckling shapes of simply supported ring- and stringer 
stiffened isotropic cylinders with four rings. 
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Figure 10.- 


(a) Internal stringers. (b) External stringers. 

Prebuckling shapes for simply supported stringer- stiffened isotropic cylinders with eccentrically 

applied compressive loads. 


NASA-Langley, 1968 32 L- 5084 
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(a) Internal stringers. (b) External stringers. 

Figure 11.- Buckling shapes for simply supported stringer- stiffened isotropic cylinders with eccentrically applied 

compressive loads. 
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u The aeronautical and space activities of the United States shall he 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 
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